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Abstract 

The mechanism of appearance of exponentially large number of metastable states in magnetic 
phases of disordered Ising magnets with short-range random exchange is suggested. It is based on 
the assumption that transitions into inhomogeneous magnetic phases results from the condensation 
of macroscopically large number of sparse delocalized modes near the localization threshold. The 
properties of metastable states in random magnets with zero ground state magnetization (dilute an- 
tiferromagnet, binary spin glass, dilute ferromagnet with dipole interaction) has been obtained in 
framework of this mechanism using variant of mean-field approximation. The relations between the 
characteristics of slow nonequilibrium processes in magnetic phases (such as hysteresis loop form, 
thermo-remainent and isothermal remainent magnetizations, field-cooled and zero-field-cooled ther- 
modynamic quantities) and thermodynamic parameters of metastable states are established. 

1 Introduction 

The specific feature of disordered magnets is the appearance of large number of metastable states in 
magnetic phases. They cause various irreversibility phenomena in the reaction to the changes of temper- 
ature and external magnetic field, such as the dependence of thermodynamic parameters on the order 
and the rate of these changes, the appearance of hysteresis loops, depending on the field amplitude etc. 
|, |, |, I. These effects appear to be common to the all types of magnetic disorder - from dilute mag- 
nets with nonmagnetic impurities ^ to spin glasses existing in solid solutions of ferromagnets with 
antiferromagnets ^ ^ . This universality of nonequilibrium phenomena compels to suggest the existence 
of common mechanism responsible for the appearance of metastable states in random magnets. 

The evidences of the existence of metastable states, the number of which growth exponentially with 
the number of sites, has been obtained in numerical studies of several models of disordered magnets 
with short-range exchange |^, |[ ^ ||] . Theoretical description of the related inergodic phenomena is also 
possible mainly with the use of numerical methods, see e. g. [|l], ^, ^. 

The principal analytical results has been obtained in the Shcrrington-Kirkpatrick spin glass model 
with infinite-range interaction. It was established that in this model the region of inergodicity exists in 
the external field less than the Almeida-Thouless field Hat, [|o|, and in the framework of the replica 
symmetry breaking scheme by Parisi ||ll| and the notion of the hierarchy of macroscopic relaxation times 
the thermodynamic parameters of field-cooled and zero- field-cooled processes has been obtained [Q. 
But this is all what is known about the irreversible processes in this model. In spite of the wide application 
of such approach to other mean-field models, see e.g. it is not clear how all variety of inergodic 

phenomena resulting from transitions between metastable states could be described by the methods of 
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Meanwhile, theoretical description of irreversible processes could be made quite simple and apparent if 
thermodynamic properties of metastable states, their regions of stability and points of possible transitions 
between them would be known at H < Hat- Principally, in the mean- field models this data could 
be obtained from (nonaveraged)equations for local magnetic moments. The most popular example of 
application of this method in the theory of disordered magnets is the TAP-equations for local magnetic 
moments in the Sherrington-Kirkpatrick spin glass model [p^ . It was established that the number of their 
solutions is exponentially large |16[, yet their explicit form has not been determined by the analytical 



methods. Now it is still not clear how many of these solutions correspond to the stable states and how 
the divergent barriers between them emerge. 

Meanwhile, the results of the investigation of TAP-equations in case of large but finite interaction 
range |p7| cast some doubt on the possibility to use the infinite-range models as a valid approximation for 
description of real disordered magnets. The authors of first noticed the consequences of the qualitative 
difference of the spectrum of the gaussian orthogonal ensemble of random exchange in Sherrington- 
Kirkpatrick model with all eigenvectors delocalized and that of short-range random exchange having 
localized eigenvectors near its edges . It was shown in that in 3d spin glass with sufficiently large 
but finite interaction range, the condensation of localized modes could drastically change the nature of 
phase transition preventing the subsequent condensation of delocalized modes if localization length index 
is greater than 2/3. This case with the absence of metastable states seems to be not very common in 
real disordered magnets but this result demonstrates the necessity to take into account the qualitative 
features of spectral characteristics of short-range random exchange. 

The studies of the spectral structure of various short-range random matrix ensembles have shown 
that all they have similar features having localized eigenvectors near the edges and fractal structure of 
eigenvectors at the localization threshold |2^. Thus one can suppose that universality of the 

inergodic phenomena in the different types of disordered magnets results from the similarity of their 
spectral properties. The common mechanism of the appearance of large number of metastable states in 
such systems could be an almost simultaneous condensation of macroscopically large number of almost 
nonoverlapping fractal modes near the localization threshold. 

Here we must note that the finiteness of the interaction range does not forbid the quantitative descrip- 
tion of such mechanism in the framework of the mean-field approach. Indeed, the macroscopic number of 
condensing modes (order parameter components) greatly reduce the number of noncondensing ones (or- 
der parameter fluctuations). So we may expect that Ginzburg parameter will be essentially smaller than 
that of the homogeneous magnet of the same dimension for all types of short-range random exchange. 

On the basis of such suppositions, we attempt here to develop the phenomenological mean-field theory 
of disordered Ising magnets with the zero ground state magnetization, using some heuristic suppositions 
about the structure of the fractal eigenvectors in random exchange ensembles. In this approach the 
appearance of numerous metastable states could be naturally explained and, using some simplifying 
assumptions, it is possible to describe analytically the thermodynamic properties of these states near 
transition point. In their turn, these results allow to determine parameters of all slow irreversible pro- 
cesses and, particularly, to get the analytical expressions, describing the hysteresis loop form at arbitrary 
amplitude of external AC field as well as field and temperature dependencies of remanent magnetizations 
in the considered types of disordered magnets. 



2 Mean-field approximation for disordered Ising magnets 

The Hamiltonian of disordered Ising magnet is 

H = -IT^J.^^^Sj (1) 
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Here is the random exchange matrix. Si — ±1. In the majority of realistic models, Jij are nonzero 
only for the sites within several coordination spheres and 

I I ^ ^max 

Here we consider just such models. 

The mean-field approximation for the Hamiltonian in Eq. |l| consists in the substitution of matrix Jij 
by the projector on its largest eigenvalue 

N„ 

J^.^jJ^cfC} (2) 

a=l 

Here cf are the normalized eigenvectors of Jij corresponding to the largest (degenerate in case of multi- 
sublattice antifcrromagnet) eigenvalue J, a = 1,..., No- 
Then the Hamiltonian becomes the function of (multicomponent) order parameter 

N 



and one should only calculate the entropy 



S (77") = lnrr(5jvi/2^o_5- ^-S; 



(3) 



to get the inequilibrium thermodynamic potential 



Minimization of F {r]°') would give the equilibrium values of thermodynamic parameters, corresponding 
to the lowest minimum, and those of metastable states, corresponding to less deep minima. In particular, 
one can get the average spin values 



In some cases this approximation may give quantitative description of thermodynamics of the second 
order transitions in homogeneous magnets (except for the fluctuation region in the immediate vicinity of 
transition point). So it is rather natural to use the mean- field approximation for the disordered magnets. 
In this case it seems sufficient to average the results obtained in the framework of the described above 
scheme over random J and cf. Then the thermodynamics of disordered magnets would not qualitatively 
differ from that of ideal magnets and, in particular, there would not appear numerous metastable states 
in the inhomogeneous magnetic phases. 

The most probable reason for the emergence of exponentially large number of metastable states lies 
in the specific structure of spectrum and eigenvectors of random matrix J^. Indeed, the eigenvectors 
of Jij having the described above properties are localized near the upper (and lower) boundary of the 
spectrum [Q. So to describe the transition in random magnet one should not take in the Eq. |^ the 
largest eigenvalue but rather lower one J at the localization threshold, i. e. the largest eigenvalue from 
those having delocalized eigenvectors. The reason is that the macroscopic transition could take place only 
as a result of condensation of delocalized mode, while the preceding condensation of local modes with 
larger eigenvalues (transition temperatures) results in the specific transition into Griffiths' phase, which 
is not accompanied by noticeable anomalies of thermodynamic parameters |2^, Yet we must note 
that according to in case of sufficiently large interaction range the condensation of local modes could 
stabilize the delocalized modes and macroscopic transition would not take place. So here we assume the 
range of interaction to be sufficiently small for macroscopic transition to occur. 

Meanwhile, it seems rather probable that in the majority of random realizations the condensation of 
just one delocalized mode is not sufficient for the stabilization of new magnetic phase. For the dimension 
d > 2, delocalized eigenvectors near the localization threshold has rather sparse (fractal) structure con- 
sisting of rare localization regions connected just by the branching chains. In other words, there is the 
set of sites with the structure of percolation cluster where N{cf)'^ » 1, while on the other sites 
N{cf Y << 1- III this respect the modes which are close to the localization threshold differ essentially 
from those in the interior of the spectrum and from the modes of translationally invariant which have 
N{cfY ~ 1 at almost all sites. 

The evidences of such fractal structure of these nearly localized modes were obtained in the numerical 
studies of various ensembles of short-range random matrices, see jl^, ^ and references therein. According 
to Eq. ^ the condensation of one such mode results in the appearance of sufficiently large average spins 
only on a sparse fractal structure, which would not suffice to stabilize the modes with lower eigenvalues, 
being localized, in general, on the other fractal sets of sites. To be more precise, the condensation of 
the nearest to the localization threshold mode can stabilize only those modes with < J, which 
overlap essentially with it, i. e. having Ncfc^ >> 1 at almost all sites where N{cf )'^ >> 1. So after 
the condensation of the first sparse mode, the second mode having almost zero overlap with the first one 
will condense at lower temperature. Further decreasing of temperature will result in the condensation 
of third sparse mode which does not essentially overlap with the first and second ones and so on. This 
subsequent condensation of almost nonoverlapping modes with lower eigenvalues will take place until 
sufficiently large average spins appear at almost all sites. In the intervals between the eigenvalues of 
condensing modes there can exist, in general, an arbitrary numbers of modes which do not condense due 
to the large overlap with the previously condensed ones. These modes represents the order parameter 
fluctuations and should be omitted in the mean-field approximation. 

Fractal structure of condensing modes suggests that their number diverges when N ^ 00. Indeed, if 
the sets of sites, where the modes considered are mainly localized, have the fractal dimension df < d, 
then the number of the condensing modes, Nq, is of the order N^~'^i/'^ . 




(4) 



The described mechanism of the transition into inhomogeneous magnetic phases can rather naturally 
explain the appearance of exponentially large number of metastable states. Indeed, the condensation of 
one mode in zero field gives rise to two stable states related by the global spin reversal, while each sub- 
sequent condensation multiplies this number by factor two. Thus there appears 2^° ~ exp(iV^~''-f/'' ln2) 
stable states. Spin configurations in these states will be related by the overturns of independent groups of 
spins, corresponding to fractal modes. Just this structure of the set of ground states has been revealed in 
recent numerical study of 3d Ising spin glass with ± J exchange [ p^ , which makes this mechanism rather 
probable. 

The fact that overlaps of condensing modes almost vanish (i. e. Ncfc'^ << 1 for a ^ (3 for almost 
all sites) allows to simplify essentially further considerations. Let us approximate the corresponding 
eigenvectors cf by the set of nonoverlapping (normalized) vectors ef, efef = for a ^ (3, which are 
equal to cf in the regions where they are mainly localized (N(cf)'^ >> 1) and zero outside them. Then 
on the subspace, spread by the (apparently, orthogonal) vectors ef , Jij is almost diagonal 

N„ 
a, (3=1 

Here Ja/3 is small positively defined matrix, 

\Jaf3\ « J 

Then it is easy to find the mean-field thermodynamic potential, depending on the multicomponcnt order 
parameter 

i 

iV„ = ^0(|ef|) 

i 

(9 (x) is Haviside's step function) and (quasi)local magnetizations 

ma = N-'J2S,0i\e:\) 

i 

It has the form 

^ No No 

^ = "2 II ^/N^{JSa,p-J^f3)U[3~TY,N^So.iL,m^)~NHm (5) 

a, (3=1 a = l 



5'a(/a,"la) =iV„MnTrc,(5j^l/2;^_^^„5^^^^,„^^^5^3(|^„|) (6) 

Here Tr^ denotes the sum over spin configuration of those sites where ef 7^ 0. According to the above 
considerations, Na — s- cx) in thermodynamic limit. 

No 
a=l 



and homogeneous magnetization is 



^" AT 

a=l 



Thermodynamic potential F, Eq. ||, depends on a small random matrix Jaf3 and random vectors ef . 
Their form is determined by the type of the random exchange matrix ensemble. In some cases it is 
possible to get some notion on the ef form. For example, in the spin glass with binary random exchange. 



in every bond configuration there are nonfrustrated d-dimensional clusters, i. e. the clusters which have 
unique spin configuration ai providing the energy minimum, and 



Thus the delocahzed eigenvectors with largest eigenvalues could be approximately constructed via con- 
nection of some nonfrustrated d-dimcnsional clusters by the branching chains without loops which are 
also nonfrustrated at all bond configurations 124]. So ef can be approximately represented as 



-1/2, 



(7) 



where af are the spin configuration constructed as described on the nonfrustrated fractal sets of sites. In 
dilute magnets with the concentration of magnetic atoms above the percolation threshold, the vectors ef 
can be also represented in the form Eq. ^ as one can connect by chains the d-dimensional ferromagnetic 
(antiferromagnetic) clusters belonging to the percolation cluster. So cr" = 1 in dilute ferromagnet and 
af = (— 1)*"^' in dilute antiferromagnet. 

Yet we must note that in some specific random realization of matrix Jy its eigenvectors near the 
localization threshold and ef could essentially differ from that of Eq. ^ Nevertheless we will suppose 
that this expression approximates ef reasonably well in the majority of random realizations so it can be 
used for the estimates of the ensemble average of the sums 



(8) 



As we show below, thermodynamics near the transition point does not depend on the detailed form of 
ef being determined by the several constants u„q,, Eq. ||. 

The advantage of F representations in the form of Eq. || is the additivity of the entropy. But 
with arbitrary random matrix J^p this expression is still difficult to analyze. We can simplify it using 
rather apparent consideration that eigenvalues of this matrix, lying in an interval between zero and some 
Jq << J, must condense near zero as the farther from localization threshold the more rare are the modes 
which do not overlap with the preceding ones. Then, as the average interval between eigenvalues is of the 
order the same order could have the smallest eigenvalues of Jaf3- So we can approximate Jaf3 by 

the projector on some (random) vector which properties are determined by the type of Jij ensemble, 



Jaf3 = JoTarp, 



No 



(9) 



Here we must note that the assumption that only one of eigenvalues of Jaf3 is finite, while all others are of 
the order 1/N, is rather rough. It results in merging of the condensation points of all modes except one 
and makes the transition to be more sharp, while actually some modes will condense somewhere between 
T = J and T = J — Jq. Yet the approximation in Eq. |9| allows to obtain analytical results which agree 
qualitatively with experiments, so it could be a starting point for more precise theory, accounting for 
distribution of the condensation temperatures of fractal modes. 

Further we will show that the form of could be determined from the fact that matrix Jaf3 defines 
the type of the ensemble's ground state. 

Then Eq. || becomes 

N„ J No No 

F = -^J2 + y/Kr^Lf -TY^NaSa (/a, m„) ~ NHm (10) 

Q — 1 a— 1 a — 1 

Partial entropies Sa {la,i^a) in Eq. ^ can be represented as 



Sa {^ai ^a) = hi 2 — max 



The values ipa , ipa corresponding to the maximum are determined by the equations 



(11) 



Differentiating potential in Eq. over and rUa we get the equations of state 

/3 



(12) 



The stable solutions of Eqs. [TTI - corresponding to the minima of F must have the positively defined 
Hessian 



Gap — Sap { T 



l-II(e?)'<^'(¥'a+V^ae?yA^) -j\+Jorarp (13) 



For = , r = we have from Eqs. |T| - |T| 

ma = N"^ ^ sign {e"la) 

i 

\la\^N-'/^Y.\<\ 

i 

Thus there are 2^" stable solutions of Eqs. |l^ - |l^ in this case which differ by the la signs . If the 
ensemble of random Jij has ground states with to = in almost all realizations then we may provide 

No 

the minimal energy for the states with m = putting ^ \/N^rala — cm. There is the unique function 

a=l 

rQ.(ef ) obeying this condition for arbitrary la signs, 

i 

Here c' is a normalization constant. 

Further we consider only disordered magnets with zero magnetization in ground states such as dilute 
antiferromagnets, spin glasses and dilute ferromagnets with dipole interaction According to the 

above considerations on the form of e" in dilute magnets and binary spin glasses (see Eq. 0), Eq. |lj can 
be represented as 

ra^U^a{Na/Y,Npu\pf'^ (15) 


Thus, in the approach outlined, the study of metastable states in the mentioned above magnets consists 
in the finding of stable solutions of Eqs. |l^ - with in the form of Eq. 15, Spin configurations 
corresponding to the obtained la^ ma can be found from the expression 



{S^)T -T.^a"'et^^^-^^+Y,N-'6{K\) (16) 

a ^ '^la a "l" 

The averaging over disorder reduces to the averaging of solutions over random ef and Jo << J- Note 
that the localization threshold J is not random quantity being the characteristics of the whole ensemble 
of random Jij . 

Smallness of Jq > means that corresponding distribution function must have sufficiently narrow 
bounded support, i. e. the possible Jq values must be smaller than some J > obeying the condition 
J « J. Contrary to the case of the sums of macroscopic numbers of variables Una in Eq. H, there are no 
reason to suppose the fiuctuations of Jo to be self-averaging, that is, that ( Jg ) — » (Jo) when N — > oo. 
So the thermodynamic parameters of metastable states in the inhomogeneous magnetic phases will not 
be the self-averaging quantities being determined by the different Jq values in different samples. Note 
that the absence of the self-averaging property of the stable thermodynamic parameters has been also 
observed in the numerical studies of disordered magnets 1^, . 



3 Thermodynamics near transition 

In the absence of external field, the equations of state, Eqs. ^ - have unique paramagnetic solution 
at T > J and a number of stable solutions appears at T < J. Thus at T = J, iJ = a transition from 
paramagnetic phase into inhomogeneous magnetic one takes place. Let us consider thermodynamics in 
the vicinity of this transition which is defined by the condition 



la, ma « 1 



(17) 



In this case the equations for magnetizations of condensing modes, nia, follow from Eqs. |Tl| - [l^ , |l5| 

Trria + TQuj^mN/ ^ iV^u?/3 + U4aml^/3ul^ = uf^H/ J (18) 

and la can be expressed via nia, 

UiJa = TOa + (uL - ^)H/J + ("Sa - "10^40) »7l^/3Ui„ (19) 

Here r = 1 — J/T, tq = Jq/T. Hessian, Eq. |l^, has the form 

From Eqs. |l^ - |l| it follows that H « J,t « 1, tq << 1. 

It is natural to suppose the sums of macroscopic numbers of variables Una in Eq- 
quantities, so we can substitute them by their average values. Assuming that Eqs, 
disorder realizations, we get 

U4a — 1, U3a = Mia, 

Let us also suppose that ufa do not depend on the a, 

_ No 

<-zi?^iVo-i^< (20) 

It seems that such approximation could not qualitatively change results. Giving it up one will just add 
some fluctuations to final expressions. Yet it allows to simplify essentially Eqs. |l^ making possible their 

analytical study. 

The constant can be estimated using Eq. |^. From Eqs. ^,|| it follows 



to be self-averaging 
hold for the most 



4 = ^0-^ E 



where i/^ and are relative parts of positive and negative values of ef , so < 1 . Then = 1 in 
dilute dipole ferromagnets. In dilute antiferromagnet the difference Va ~ can be nonzero only due to 
uncompensated spins on the surface of d-dimensional antiferromagnetic clusters on which ef are mostly 
localized. Hence Va — i^a of the order of the surface to volume ratio of d-dimensional antiferromagnetic 
clusters, so ~ D~^, where D is the average diameter (in terms of lattice spacing) of these clusters. 
Evidently, D is a function of the concentration of antiferromagnetic atoms, which goes to infinity when 
concentration tends to 1. 

In the binary spin glass, uf depends on the concentration of ferromagnetic bonds p. In this case 
Ui = l for p > I — Pc and = for p < pc, Pc being the bond percolation threshold on the lattice of 
magnetic atoms. At < P < 1 ^ Pa u\ dependence on p can be qualitatively described as 

— p-pc + D-^{l-p-pc) 
= 

Here we also substitute Na by their average value 

Na = N/No. 

Then, introducing instead of nia reduced magnetizations , 



and dimensionless field 



Eqs. EG can be represented as 



h = \[^^H/T 

Tfla+ToH + nl/3^h (21) 



where p. = ^SaMa- 



At r > Eq. ^ have unique paramagnetic solution with equal — A*: which we denote as fiQ- It 
obeys the equation 

{T + To)no+f^l/3^h (22) 
When T < Eq. ^ have 2^" — 2 stable inhomogeneous solutions besides which can be represented as 

fJ-a = V—T (sin ip + VSaa cos ip^ . 
Here era = ±1 and fir, tq, h, A) is the solution of the equation 

3to (VsA cos If + sin tp^ — 2r sin 3<y9 = 3/i |t| , (23) 



(24) 



The parameter A, varying in the interval (—1, 1), defines the degree of the inhomogeneity of a metastable 
state and A = ±1 correspond to the paramagnetic state with /i^ = M- states with equal A have also 
the equal magnetizations 



/i = ^/—T (sin if + y/SA cos ip^ , 
and Edwards- Anderson order parameter 

q^N-^Yl (^')t - ~ II - = 3t(A2 - 1) cos2 + (1 - ^) 

i a 

as well as thermodynamic potential 

4F/TN = (ro + itfr)^^ + rq - Shfj. - 4 In 2 
These states are stable at r + > 0, which is equivalent to the inequality 

\ip\ < 7r/6 



(25) 



(26) 



(27) 



(28) 



As To > 0, the left side of Eq. 23 is a monotonously growing function of ip for \ip\ < 7r/6. Hence, there 
is only one stable solution for ip at a, given A, which exists in the interval of the fields 

h_<h<h+, h±^ V3hATA/2±hc 



Hat = V-3Tro 



(29) 



he = \/^(to/2 - 2r/3) (30) 
In this interval the solution of the Eq. |2^ can be approximated by the quadratic function of the field 

4(2 - V3)AhAT{h+ - h){h - h-) 



if : 



2h - V^HatA 



^hl - (2 - %/3)2A2/i2 



AT 



(31) 



which gives the exact values ip{h±) = ±7r/6 and p{AhAT) = 0. 

From the stability condition, Eq. ^ and Eqs. ^ ^ it follows that metastable states are stable in 
the region 

9(ToAi - hf < -4t3 

which is the band on the ^, h plane. The magnetization is a monotonously growing function of h and A 
inside this band so the field dependencies of magnetization can be represented as a set of uncrossed lines 
bounded from above and below by the fJ-o{h) curve as shown in Fig. 1. 

Apparently, the region on this figure, where metastable states exist, defines the form of hysteresis 
loop, which appears as a reaction on a slow AC field with amplitude greater than 



he — he + V^hAx/'^ 



(32) 



Figure 1: Field dependencies of magnetization in metastable states (dashed lines) and stable states (solid 
lines) near transition, (a) - < — 2t < 3to, (b) - 3ro < — 2t. 



There are certain temperature variations in the loop form as at — 2t < 3to only part of metastable states 
are stable a.t h — and loop is rather slim, sec Fig. 1(a), while at — 2t > 3to all metastable states are 
stable at h = and loop became more thick, Fig. 1(b). Note also that when amplitude of AC field is less 
than he the form of hysteresis loop is defined by the field dependencies of magnetizations in corresponding 
metastable states. 

Let us also present the expressions for (dimensionless) magnetic susceptibility x = dfi/dh, 



X ^ = To — 2rcos3(p/(cos(p — V^Asintp), 



entropy S and heat capacity C, 



C= 1 



■3 ro(l-A2) 



X 



2 1 - Vi^tgip 



To 



The above expressions allow to get some notion about the field and temperature dependencies of these 
quantities. Thus at the boundaries of stability region, h = h±, q and x^^ has the lowest values 



q = 9t (A2 - 1) /4 - (1 - ufj T (3A ± 1)' /4, X = 1/to, 
while magnetization, entropy and heat capacity are 

= V^(3A± l)/2, 5 = ln2 + r(5±3A)/4, C = 3 (1 ± A) /2 - t/to 



(33) 



(34) 



When \h\ goes to h^, Eq. the more homogeneous states with A ^ ±1 stay stable and their magne- 
tization tends to fj,o (±/ie) = ±2^/—t. However the limiting values of magnetic susceptibility and heat 
capacity differ from those in paramagnetic state: 

X(7^ =T + To + /io, Co = (t + To + /io) . 
In the middle of the stability band {at (p — Q or h — Ah at) we get: 



H = AT^, q = -3t ( 1 - ujA^ 



X 



ro-2r, S' = ln2 + 3r/2, C 



1- A^ 



To 



T0-2T 



In this case with the diminishing of inhomogeneity when A — > ±1 or ft, — > ±hAT, Xj ^^'^ ^ tend to 
the corresponding values of the paramagnetic phase. 

The Almeida- Thouless field /iat, Eq. 29, determines (to the order A^q"^) the point of the transition 
into the paramagnetic phase. To show this let us find the values A^q corresponding to the states with 

we get 



the lowest potential. Differentiating F , Eq. |27|, over A and using Eqs. |23|, |25| and Eq, 



dF . 
— — = —A^T sines cos 09, 
dA ^ ^ 



5A2 



0=0 > 



Thus the lowest value of potential have the states with (/3 = at given r and h . One can see that Eq. 
has solution with ip — when A — h/hAT which is possible at /i^ < hj^j,. When h'^ > h'^j, F (A) has no 
minima inside the region A^ < 1 in which it is defined and the minimal values occur at its boundaries 
for Aeq = sign{h). So the transition into paramagnetic state takes place at /i = ±hAT- 

As A is a rational number of the form 2n/No — 1 (cf. Eq. ^J) it can not be exactly equal to h/hAT 
at all h^ < h\rp. Hence A^q is defined so that |A — h/hAT\ is minimal and can be represented as 



JVo-l 



n=l 







AT 



h 2n 

hAT No 

Thus at h^ < h\rp series of transitions between inhomogeneous magnetic states takes place at fields 

hn = hAT ( — TT— - 1 
The value of ifeq at /i^ < /i^j. corresponding to Aeq is, see Eq. GS 



/o No-1 



n=l 

Inserting Aeq and ifeq in the Eqs. ^we get the equilibrium values fieq and qeq for N — s- oo: 

To ' 

qeq = -3t (1 - h^hlj.) e {h\T - h^) + (l - ^) Me, 
Differentiating rrieq over h we get the equilibrium susceptibility 

Xeq = To^O {hlr ~ h^) + {r + To+ Mo)"' {h^ - h^^) . 

The equilibrium entropy can be obtained by the differentiation of the equilibrium potential which to the 
order is 

Feq ^F{A^ h/hAT) - TSconf 

where configurational entropy Sconf is determined by the logarithm of the number of states with the 
same potential F, 

Scojif is of the order Nq/N and can be neglected. Hence at A^o ^ oo 



Seq = \ii2 + ^0 {h\^ - h') - f^e {e - h\^) 



For the equilibrium heat capacity we get 

= i^^AT - — 2 ^ {^'^ ~ ^AT) 

^ T + To + Mo 



Let us note that /Zgg, qeq and Seq are continuous at h? = /i^t^, while Xe, and Ceq undergoes jumps at the 
transition into paramagnetic phase. 

We must also note that the average equilibrium parameters are generally unobservable quantities due 
to the macroscopic free energy barriers between metastable states. Probably, the experimental values, 
which are rather close to them, are obtained after cooling in small external fields down to T just below 
J (field-cooled (EC) regime) ^ when barriers between metastable states are relatively small and 
system could relax into the lowest (or close to it) state at a sufficiently slow cooling. In zero field cooled 
(ZEC) regime when field is applied after cooling below T = J in zero field, the observed thermodynamic 
parameters would differ from equilibrium ones as the system would at first be trapped in the state with 
A = and will stay in it if applied field does not exceed /ic, cf. Eq. M and Eig.l. Thus at h < he the 
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Figure 2: Field (a) and temperature (b) dependencies of some thermodynamic parameters in FC and ZFC 
regimes, 1-FC susceptibility, I'-ZFC susceptibility, 2-FC magnetization, 2'-ZFC magnetization, 3-FC EA 
order parameter, 3'-ZFC EA order parameter 



ZFC parameters are those of A = metastable states. Their values can be obtained from the general 
expressions at A = and (p — irh/Ghc, see Eq. When applied field h > the system relaxes into 
the metastable state at the boundary of stability region (on the lower branch of hysteresis loop) with 

AzFC = 2{h- hc)/V3hAT. 

Inserting this Azfc in Eq. ^ and Eq. ^ (with plus sign) we get the values of thermodynamic param- 
eters the observed quantities would relax to in ZFC regime at hf. < h < h^. And when h > ZFC 
parameters correspond to those of paramagnetic state. The field and temperature dependencies of some 
thermodynamic parameters in FC and ZFC regimes are shown in Fig. 2. 

Similarly, the regions of existence of metastable states, see Fig. 1, and their parameters define the 
other quantities which are determined in the slow nonequilibrium processes, such as thermo-remainent 
magnetization, ^trm, which remains after FC process and subsequent switching off the field, and isother- 
mal remainent magnetization, ^lirm, remaining after ZFC process followed by the application for some 
time (longer than the intravalley relaxation time) an external fiel d [P, ^. Apparently, ^irm is nonzero 
only at h > he and an expression for it can be obtained from Eq. g5| at 

AiRM = min l,2hc/V3hAT,'2{h ~ hc)/V3hAT 

II, n \ \ V^T^hATAjiiM 
^IRM = (p{h = 0, AiRM) ~ ^ 

IJ-iRM can also be obtained from Eq. |25| with 

Atrm = min (l, 2hc/V3hAT, h/hArj 

VSnhArATRM 



I2hr 



The field dependencies oi ^trm and (Iirm are shown in Fig. 3(a). At ft > Htrm = min{2hc/ V3, hAx) 
IJ'TRM becomes constant, while above hiuM = min{2hc, he) ^iirm also saturates at the same value. This 
value of saturation magnetization is 

Aioo = e(3ro + 2T)2(-T)3/73ro + h3(ro + T)]^''^e{-iTQ - 2t) 



Temperature dependence of /ioo is shown in Fig. 3(b) 




Figure 3: (a) - field dependencies of ^trm (1) and hirm (2), (b) - temperature dependence of saturation 
magnetization /ioo 



4 Conclusions 

The main result of the present work consists in the qualitative, but complete description of properties of all 
metastable states in the inergodic phases of random Ising magnets with zero ground state magnetization 
and elucidation of their relations to the parameters of slow irreversible processes. The results depicted 
in Fig. 1 allow to describe every conceivable irreversible process with arbitrary sequence of field and 
temperature changes. Qualitative agreement of the obtained here parameters of some such processes 
with experiments and numerical studies jl]- justifies the approximations used in Eqs. 0, ^, |l^, 
and shows that the condensation of macroscopic number of sparse fractal modes near the localization 
threshold do can be a possible mechanism of the appearance of exponentially large number of metastable 
states in disordered Ising magnets. The argument in favor of this mechanism is also the structure of the 
set of these states related via overturns of independent spin groups, corresponding to fractal modes, as 
just the same relations between ground states in Ising spin glass with ± J exchange are revealed in 
recent numeric studies §]. 

Let us also note that obtained here results are expressed solely in terms of statistical characteristics 
of random exchange matrix. So the present approach could serve as a starting point for the developing 
of more precise quantitative theory of metastable states in disordered magnets. Such theory should be 
based on the detailed studies of the properties of random Jij eigenvectors near the localization threshold, 
which we were compelled to describe here in terms of phenomenological suppositions. One of the tasks of 
this theory could be, in particular, the test of the universality of the properties of the magnets with zero 
ground state magnetization, as the results for them differ in present approach just by the values of u\. 

This work was made under support from Russian Foundation for Basic Researches, Grant N 98-02- 
18069. 



References 

[1] K. Binder, A. P. Young, Rev. Mod. Phys., 58 (1986) 801. 

[2] R.B. Stinchcomb in Phase Transitions and Critical Phenomena, v. 7, ed. C.Domb, J.L.Lebovitz, 
Academic Press, London (1983). 

[3] C. Y. Huang, J. Magn. Magn. Mat., 51 (1985) 1. 

[4] C. Ro, G. S. Grest, CM. Soukoulis, Phys. Rev. B 31 (1985) 1682. 

[5] P. Sibani, preprint |cond-mat/9706232| (1997). 

[6] S. Bastea, P. M. Duxbury, preprint |[:ond-mat/9801108| (1998). 

[7] A. K. Hartmann, preprint |c;ond-mat/9902120| (1999). 



G. Hed, A. K. Hartmann, D. Stauffer, E. Domany, preprint [cond-mat/0007356| (2000) 



[9] A. A. Likhachev, preprint |cond-mat/0007504| (2000) 



[10] J. R. L. de Almeida, D. J. Thouless, J. Phys. A 11 (1978) 983. 

[11] G. Parisi, J. Phys. A 13 (1980) L115, 1101, 1807. 

[12] H. Sompolinsky, Phys. Rev. Lett. 47 (1981) 935. 

[13] C. De Dominicis, Y. Y. Goldschmidt, Phys. Rev. B 41 (1990) 2184 

[14] A. Crisanti, H.-J. Sommcrs, Z. Phys. B 87 (1992) 341. 

[15] D. J. Thouless, P. W. Anderson, R. G. Palmer, Phil. Mag. 35 (1977) 593. 

[16] A. J. Bray, M. A. Moore, J. Phys. C 13 (1980) L469. 

[17] L. B. Joffe, M. V. Feigelman, Zh. Eksp. Teor. Fiz. 89 (1985) 654. 

[18] I. M. Lifshitz, S. A. Gredeskul, L. A.Pastur, Introduction to the theory of disordered systems. Wiley, 
New York (1988). 

[19] M. Janssen, Int. J. Mod. Phys. B 8 (1994) 943; Phys. Rep., 295 (1998) 1. 
[20] F. Evers, A. D. Mirlin, preprint Fond-mat/0001086| (2000). 



[21] R. B. Griffiths, Phys. Rev. Lett., 23 (1969) 17. 
[22] A. J. Bray, Phys. Rev. Lett., 59 (1987) 586 . 

[23] D. Stauffer, A. Aharony, Introduction to percolation theory, Taylor&Francis, London (1992). 

[24] A. Coniglio, F. di Liberto, G. Montroy, F. Peruggi, Phys. Rev. B 44 (1981) 1665. 

[25] M. J. Stephen, A. Aharony, J. Phys. C 14 (1981) 1665. 

[26] S. Wiseman, E. Domany, preprint |cond-mat/9802095| (1998). 



